Breathers and their interaction in the massless Gross-Neveu model 



Christian Fitzner* and Michael Thies^ 
Institut fur Theoretische Physik III, Umversitat Erlangen-Niirnberg, D-91058 Erlangen, Germany 

(Dated: October 17, 2012) 

The breather is a vibrating multifermion bound state of the massless Gross-Neveu model, origi- 
nally found by Dashen, Hasslacher and Neveu in the large N limit. We exhibit the salient features 
of this state and confirm that it solves the relativistic time-dependent Hartree-Fock equations. We 
then solve the scattering problem of two breathers with arbitrary internal parameters and veloci- 
ties, generalizing an ansatz recently developed for the baryon-baryon scattering problem in the same 
"^Cg". model. The exact analytical solution is given and illustrated with a few examples. 
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O 1 I. INTRODUCTION 

o : 

The principal sources of experimental information about strong interactions are hadron spectroscopy and hadronic 
scattering processes at accelerators. Theoretically, spectroscopy is by now fairly well understood owing to lattice 
simulations of quantum chromodynamics, at least for the low-lying, most stable states. Since it is necessary to work 
in Euclidean time when computing the path integral by Monte-Carlo methods, this tool fails in the case of scattering 
problems. As a consequence, there has been little progress towards understanding scattering of composite, relativistic 
JL . bound states in the non-perturbative regime over the last decades. Early attempts to use hadron models of confined 
quarks for scattering were plagued by difficulties with covariance. It is expected that only a relativistic quantum field 
theory can account correctly for covariance. Then one is immediately faced with the lack of calculational tools at 
strong coupling in Minkowski space. This is one of the motivations for us to step back and study scattering processes 
of composite objects in an exactly solvable toy model, the Gross-Neveu (GN) model 

As is well known, the GN model is the 1+1 dimensional relativistic quantum field theory of N flavors of massless 
" Dirac fermions, interacting through a scalar-scalar contact interaction, 
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" A number of exact, analytical results about the static properties of this model have been obtained notably in the 
^—1 . 't Hooft limit (N — > oo, Ng 2 — const.), ranging from the rich hadron spectrum and hadron structure @, Q to the 
phase diagram at finite temperature and chemical potential During the past few years, time-dependent issues 
. £h ■ boosted hadrons and structure functions [j| or scattering processes involving kinks [(| 0] , kink-antikink baryons 

Q and composite multibaryon states || have been solved exactly. Although a full mathematical proof of the recent 
results in [9| is still missing, it is probably fair to say that we understand the scattering of any number of boosted, 
static hadrons of arbitrary complexity in the initial and final states. 

As pointed out in this cannot be the whole story, even in such a simple toy model as the GN model. Dashen, 
Hasslacher and Neveu (DHN) have already discovered a "breather" solution long ago i.e., a multifermion bound 
state oscillating in time in its rest frame. It generalizes the ground state baryon to a collective excited state. Breathers 
are a well-studied soliton species in the field of nonlinear science, the best known example probably being the sine- 
Gordon breather, a vibrating kink-antikink state. The theoretical interest in this particular kind of solitons stems 
from the role they play in many different areas, such as condensed matter physics, hydrodynamics and nonlinear 
optics, see e.g. |10l - fl4| and references therein. (One particularly nice application is to view the breather as a moving, 
relativistic clock for the twin paradox of special relativity [HI-) Breathers are less familiar in particle physics, as 
they have no obvious correspondence in the known particle zoo. The reason lies in their genuine classical character. 
In the GN model, they owe their existence to the large N limit, an idealization which is quite far from reality in 
particle physics. A useful analogy in strong interaction physics can nevertheless be identified, namely the collective 
vibration of a heavy nucleus. In the large N limit, hadrons may be thought of as systems made of a large number of 
constituents which can then exhibit classical behavior, much like heavy nuclei or molecules in the real world. 
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The main difference between the DHN breather and the sine-Gordon breather is the fact that the GN model is 
a fermionic theory. Here, the breather describes the mean field (Hartree-Fock potential), generated dynamically by 
quantized fermions that it drags along. We cannot start from a classical, bosonic, non-linear equation like the sine- 
Gordon equation, but have to solve a quantum mechanical self-consistency problem, including the polarization of the 
Dirac sea. Describing how one can find such solutions systematically and clarifying the role of the fermions populating 
the breather, also in a breather-breather scattering process, will be the main topics of this paper. 

In Sect, m we will review in detail the DHN breather. This particular solution of the GN model has never been 
discussed in any detail, to the best of our knowledge. In Ref. @ , it has only served to derive kink-antikink scattering 
by analytic continuation in one parameter, following a suggestion of DHN in their original work Q. The main part 
of the paper, Sect. IIII1 is then devoted to the intricate scattering process of two breathers. This covers at the same 
time breather-breather bound states as well as all problems that can be arrived at by replacing one or both of the 
breathers by a DHN baryon or a kink. This whole class of problems will be solved exactly with a generalization of 
an ansatz method developed for baryon-baryon scattering in [8|. We finish with a brief summary and an outlook, 
Sect. EH 



II. SINGLE DHN BREATHER 



Before attacking the complicated breather-breather scattering problem, we have to get familiar with the breather of 
the GN model originally found by DHN Q. The present section serves to introduce the breather, illustrate its salient 
features and show how to generalize the ansatz technique of Ref. Q to breather- type solutions. We also have to fill 
some gaps in the literature like computing the breather mass and its fermion density. From a technical point of view, 
we shall use this section to set up a convenient notation, a prerequisite for the more involved problems to follow. 



A. Reminder of the self-consistent scalar potential 

In their seminal paper where baryons of the GN model were first constructed 0, DHN also report on a time- 
dependent, semi-classical solution, the breather. Whereas the baryons were derived by means of inverse scattering 
theory, the form of the breather self-consistent scalar potential has been guessed by the authors, using the analogy with 
the well-known sine-Gordon breather. The parameters were then determined self-consistently. DHN also suggested 
that analytic continuation in one of the parameters to imaginary values should describe kink-antikink scattering, a 
suggestion which was taken up and verified in detail in Ref. Q . There, one can also find expressions for the breather 
spinors corrected for misprints in the original reference. In the present subsection, we recall the form of the scalar 
potential and fermion density in the original notation of DHN and illustrate it with a few representative examples. 

Following previous works, we identify DHN's semi-classical path integral method with the relativistic version of 
time-dependent Hartree-Fock (TDHF) in the canonical framework, i.e., the solution of the self-consistency problem 

occ 

W-S)1> a =0, S = -.9 2 E^V (2) 



The sum runs over all occupied single fermion states, including the Dirac sea. Units in which the vacuum fermion 
mass is 1 will be used throughout this work. DHN write the scalar potential as 

s _ coshKx-acosttt + b(l- K 2 /2) 
cosh Kx + a cos Qt + b 



There are two independent variables, b and e. The parameters fi, K can be expressed through e as 

7t 



n = .=£= , K = eil, (4) 



whereas a is the solution of the equation 

= b 2 K 4 + 4K 2 (1 - b 2 ) + 40 2 a 2 . (5) 

We may choose the positive square root without loss of generality, since the other sign merely corresponds to a shift 
of t by half a period. Thus S is fully specified by the parameters e > and b > y/l + e 2 , where the latter bound 
follows from ([5]). Just like the DHN baryon, the breather has two bound states. The lower one is taken to be fully 
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occupied, whereas the upper one carries vN fermions (v = 0...1 will be referred to as occupation fraction). The 
antibreather would have the upper state empty, the lower state occupied by (1 — v)N fermions, and identical S. We 
do not consider more general ways of filling the bound states which also exist, just like for baryons. DHN find the 
following self-consistency condition relating i/, b and e, 



6 = (1 _„)_^ + Z . ( 6 ) 

v 'l - (2/7r)arctane v ; 

It is worth mentioning that the breather contains the DHN baryon as a special case. For a = or, equivalently, 
b = y/l + e 2 (its minimal value), S becomes static, 

2e 2 

5=1- ; (7) 

TIT? (cosh -^m^ + VTT^ 

Comparison with the original parametrization of the DHN baryon Q , 

5 = 1 + y(tanh£_ — tanh£+), £± = yx ± — artanhy, (8) 



shows perfect agreement for the choice 



V = , (9) 



Likewise, the self-consistency condition (J6)) reduces to 

2 . . 

1 - ^stat = 1 arctane, (10) 

in agreement with the baryon case. Nothing like this could happen in the case of the sine-Gordon breather. Here, 
there is simply no static kink-antikink bound state to which the breather could possibly be reduced. This difference 
is due to the valence fermions in the GN case which overcome the kink-antikink repulsion and lead to bound, static 
kink-antikink states. 

The physical meaning of the parameters e, b can be exhibited as follows. Choose a value of e > 0. For the minimal 
allowed value of b (b = s/1 + e 2 ), the breather becomes static and e determines the size and shape of this baryon as 
well as its occupation. If we now increase b, the frequency of the breather does not change (being solely determined 
by e), but the amplitude of the oscillation increases. At the same time, the occupation fraction v decreases, 

1_ " 6 :=A>1. (11) 



1 - Vstat VH 

Let us try to give an overview of how the breather behaves in space and time. First consider the allowed range of 
parameters, using (e, A) rather than (e, b). As illustrated in Fig. [TJ the physical region in the (e, A) plane is restricted 
to 

1 < A < {^-- arctane^ . (12) 

These limits are shown together with curves of constant v in Fig. [T] Along the upper boundary, fermion number 
vanishes [y — 0). The lower boundary (A = 1) corresponds to the baryon where e and v are related according 
to Eq. (jlpp . Thus for a given occupation fraction v 1 a one-parameter family of breathers exists. By contrast, the 
sine- Gordon breather has only a single parameter governing both its frequency and amplitude and carries no fermions. 

We illustrate the dynamics of the breather for two values of e. The shape of the baryon potential evolves from a 
shallow, attractive well for small e to a widely separated kink-antikink pair at large e. The breather potential at t = 
is qualitatively similar and oscillates monotonically (but anharmonically) between two limiting curves, with a period 
of T = 2tt/Q. The fermion density of the breather is computed by adding up the contributions from the Dirac sea 
and the two bound states. If the lower bound state is completely filled, its contribution is cancelled by the density 
induced in the Dirac sea and the full density is simply given by the upper bound state, just like for the baryon. The 
total fermion number is Nv. 

Let us first illustrate the scalar potential and the fermion density for the moderate value e = 2 and two different 
values of A. For A = 1.1, v = 0.6753, S oscillates, staying always below the vacuum value S = 1, see Fig. [5] There, 
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one also sees that the density oscillates between two peaks and a single peak. If we increase the value of the breather 
amplitude by choosing A = 3.38, v — 0.002335 and the same e, S overshoots near the center of the breather, reaching 
the value 1.39, see Fig- El Between the two examples shown in Figs. [5] and [3] there is a "critical value" of A, b = 1 + e 2 
or A = yl+?i where the potential oscillates between a baryon-like shape and the constant vacuum value (S = 1). 
For this particular value of b, S becomes 

cosh Kx — e 2 cos Qt + 1 — e 2 
crit = coshKx + e 2 cosm+l + e 2 ' ( ' 

Indeed, after 1/2 period, cos Sit = — 1 and S^it = 1 for all x. This special case is noteworthy for yet another reason, 
namely its close relationship to the sine-Gordon breather. The sine-Gordon breather is a solution of the sine-Gordon 
equation 

9 M 9^ + sm<^ = (14) 

given by 



(\J\ — lJ 1 COS Celt \ , . 
1^=^^ ■ (I 5 ) 
lj cosh yl - uj 2 x I 
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FIG. 3: Like Fig. [21 but for A = 3.38 corresponding to v — 0.002335. The potential S exceeds the value 1 near the midpoint 
by 40%. Note the different scales in the density plots of Figs. [2] and [3] 



One can check that 

Sent = cos ^, (16) 

provided one identifies u) with 0/2= + e 2 . The motivation behind the particular nonlinear transformation (|16j) 

is the fact that it also maps the sine-Gordon kink, <f> = 4arctane :E , onto the Gross- Neveu antikink, S — — tanhx. We 
do not know whether this is a mere coincidence or whether there is a deeper reason behind this mapping. 

We now turn to a large value of e (e = 700) where the breather has a more pronounced kink-antikink shape. Fig. 0] 
shows plots of S and the fermion density at A = 10, v = 0.9909, where the system exhibits oscillation of a kink against 
an antikink without overshooting. If we go beyond the critical value A = yl + e 2 , S again exceeds the value of 1 




FIG. 4: Scalar potential S (left plot) and fermion density p (right plot) during one period of the breather. Parameters: 
e — 700, A = 10 corresponding to period T = 2199 and occupation fraction v = 0.9909. Kink and antikink stay well separated 
during the whole period. 

at the midpoint. However the time scale for this overshooting is now much shorter than the period of the breather, 
indicating that the motion is strongly anharmonic. For this reason it is difficult to present a plot similar to Fig. [3] 
here. Instead, we show the very short time interval where S exceeds 1 in "slow motion" , see Fig. [3] 

B. Light cone variables and rational parameters 

The DHN breather has 2 real parameters, e and b, governing its size, frequency and amplitude in the rest frame. 
In addition, one can perform arbitrary Poincare transformations, yielding 3 additional parameters (velocity, shift in 
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x and t). Breather-breather scattering then will depend on 10 parameters, as compared to 6 parameters for baryon- 
baryon scattering. In view of this complication and according to our prior experience with baryon scattering [a, 
it is crucial to simplify the kinematics as much as possible by going to light cone variables. Moreover, it is highly 
advisable to introduce parameters such that the final results are rational functions, rather than algebraic or tran- 
scendental functions. This is a prerequisite for being able to reproduce lengthy computer algebra (CA) computations 
independently and simplify them to a unique form. In this subsection, we will take over the notation from previous 
papers whenever possible and extend the rational parametrizations from the baryon to the breather. 
Our convention for light cone variables is 

z — x — t, z^x + t, do = B — d, d\=d + d. (17) 

We use the light cone spectral parameter 

rather than momentum k and energy w. The boost parameter rj, velocity v and rapidity £ are related by 



T > = #=\lT ± T.> U = V^- ( 19 ) 



' 1 + v _ rf - 1 
1 — v ' rj 2 + 1 

Under a boost, the light cone coordinates and spectral parameter are simply rescaled, 

z^rjz, z^rrf 1 ^, C — > VC- (20) 
In these variables the Lorentz invariant argument of a plane wave reads 

kx - cut = - \(z- - ] . (21) 



Together with the chiral choice of Dirac matrices (diagonal 75), 

7° = <ri. 7 1= ^ 2) 75= 7V = -<t 3 , (22) 
the Dirac equation entering the TDHF approach simplifies to 

2idip 2 = Sipi, 2i&ipi = -Sip 2 , (23) 

with chirality = ipL,ip2 = i>R- In the case of the DHN baryon, the following reparametrization of y has proven to 
be useful to avoid the appearance of square roots, 

. wi=\l-y? = i— — . 



w = ^ = ^-yi=-i^> z ^m- (24) 
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Introducing the phase tpi, 

Z\ — e lipi , yi — simpi, wi=—cos(pi, (25) 

we see that one should restrict Z\ to the 2nd quadrant (7r/2 < <p\ < 7r) in order to parametrize the relevant region 
Vi,wi > 0. These are all the ingredients without which the solution of baryon scattering problems would have been 
almost hopelessly complicated. In the same vein, let us reparametrize the breather specific quantities (K, fl, e) and 
(a, b, A). For K, f2, e, we choose again a phase Z £ U(l) via 

r, Z 2 -l n Z 2 + l K [Z 2 -l\ 

Upon setting Z = e lip , we find 



K = 2sixi(p, fi = — 2 cost/?, e = — tamp, \/l + e 2 = , (27) 

COS</? 

so that if , f2, e > corresponds to ip £ [tt/2,tt]. The remaining set of parameters (a, 6, A) require yet another phase 
variable Q £ U(l) in the 2nd quadrant. We set 



Z 2 + 1 Q 2 + 1 ' Z 2 + 1 

so that Eq. ([5]) is satisfied identically. Parametrizing Q = e 1 ^ , one finds 

a = tant^tan?/;, b— -, A = -, \/ X 2 — 1 = — tan tp, (29) 

cos (f cos "0 cos 

hence a, 6, A > is consistent with ip £ [tt/2, it]. Thus the two new parameters (Z, Q) replacing (e, b) both live on the 
quarter of the unit circle in the 2nd quadrant. The baryon limit of the breather (6 = vl + e 2 ) amounts to setting 
Q = — 1. In this limit, Z has the same meaning as in our previous works on baryon scattering [1, [§]. The value 
b = 1 + e 2 where the DHN breather can be mapped onto the sine-Gordon breather corresponds to the choice Q = Z. 



C. Breather mass 

The mass of the breather has apparently not yet been calculated in the literature. Here, we closely follow Ref. @ 
where the original computation of the baryon mass was simplified. As explained there, it is advantageous to 
compute the total energy by splitting the energy density into local (vacuum subtracted) and constant pieces. The 
local part (which gives a finite contribution when integrated over dx) consists of 3 contributions: The kinetic energy 
from continuum states, subtracting the asymptotic value 

- N L d *C ^ H*** + ■ (30) 



the subtracted potential energy 



e£H = -JL:I d x (S 2 -l), (31) 



2Ng 2 



and the kinetic energy from the discrete states 



-iN / dx (V ^Sa^W + W (2)t 75<9zV (2) ) ■ ( 32 ) 

J — oo 

Using the vacuum gap equation 



(33) 
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(2) 

and performing all integrations analytically, one finds as usual that E^ oc cancels exactly the logarithmic divergence 
in E^l . The finite part of E^l combines with E^l to a complicated expression multiplied by a factor 

2Q 4tQ\nZ 

which vanishes owing to the self-consistency condition. Hence the local energy part of the energy vanishes. The total 
mass of the breather is then completely determined by the fermion phase shifts 2, 5], which are identical to those of 
a baryon with the same value of Z or e. We conclude that the mass of the breather with parameters (b, e) is the same 
as the baryon mass with parameter e, 

2eN 

M breather (e, b) = — j= = . (35) 



ttVTT 



Thus the breathers shown in Figs. [5] and [3] have a common mass (M/N — 0.56941003), as do the breathers shown in 
Figs. |4]and[5] (M/N = 0.63661912). Independence of the mass on the breather amplitude is counter-intuitive at first 
sight. However, we have to remember that the baryon number depends on Q and decreases with growing amplitude, 
so that there are competing effects. If we compare the mass of the breather with the mass of a baryon with the same 
fermion number, using Eq. dill) , we indeed find that the breather gets heavier with increasing amplitude, 

2N (ttv\ 2N /tt(A-1 + i/)\ . , 

Mbaryon = — Sin ^ — J , A/brcathcr = — Sin I — I . (36) 

At fixed v, the ratio Mb rea ther /-^baryon increases monotonically with A. Clearly, the breather mass always stays below 
the kink-antikink threshold 2N/tt, a precondition for its stability. 



D. Exponential ansatz method for single breather 

Scattering problems involving baryons and multi-baryon states have been solved recently by means of a joint ansatz 
for the scalar TDHF potential and the Dirac spinors Assuming that the self-consistent potentials are transparent 

and working with rational functions of certain basis exponentials, it was possible to solve the Dirac equation purely 
algebraically. This method has the potential to handle the breather-breather scattering problem as well. In the 
present section, we cast the results for the single breather into a form well suited for such a scattering calculation. To 
this end, we have to express the single breather potential and spinors in terms of basis exponentials. To account for 
the new fact that both hyperbolic and trigonometric functions appear in S, Eq. ([3]), we evidently must complexify 
the basis exponentials. 

This section is also a preparation for tackling the breather-breather scattering problem in Sect. IIIII There, we 
obviously will have to deal with 2 different breathers, moving at different velocities. To keep the notation consistent, 
it is therefore advisable to introduce breather labels i = 1,2 right away, even if it renders the notation in the present 
subsection more cumbersome. To distinguish the breather labels from many other subscripts or superscripts, we will 
use superscripts (i) with parentheses. 

So far, we have only considered the breather at rest. Boosting the breather to different velocities is straightforward 
in light cone coordinates. We will indicate the necessary modifications introduced by the boost as we go along. 

We introduce as basis exponentials for breather i in the notation of subsection III Bl 

^ (i) =ex P (g-|f), VP = VP\ (37) 

The boost parameter rji depends on the breather velocity Vi, see Eq. (|19|). The breather potential (|3J) can be written 
as the following rational function of Vy , , 



= 1 + a»(V«) 2 + agvfVW + a«(vf ) 2 + a« ,{V^f{V^f, 
2>W = 1 + ^(Vf ) 2 + ft&V^VW + &g(vf ) 2 + b^Uv^nV^f, (38) 



9 



with 

_(<) _ „(i) _ _ _i.(0 _ Z i ~ 1 gj ~ 1 

Oil " «22 - "I! " »22 - Z 2 + 1 g 2 + r 

(t) _ ^ 4 + 1 2Z a 2Q t , (i) _ 2Zj 2Qj 

° 12 ~ Z 2 £? + lQf + l' 12 ~ Z? + 1Q? + 1' 

a 1122 = ^1122 = 1- (39) 

The ansatz for the continuum spinors (assuming a transparent potential) follows exactly the strategy used earlier for 
baryons 8], 



with 



A/f = i + c «(y«) 2 + 4V«v« + c«(y«) 2 + c« 2 (^«) 2 (k«) 2 , 

A/f = i + 4? (y«) 2 + dgv^vf + 4 } (vf ) 2 + $} m (y®f{v®)\ (4i) 

The free spinor has been pulled out so that all polynomials start with a 1. When inserted into the Dirac equation, 
this ansatz leads to an overdetermined algebraic problem. We find a unique solution with coefficients 

.(») _ _M _ + 

(*) _ _j(») _ C^gi - 1 (0 

C22 ~ 22 ~ CviZi + i 11 ' 

(i) _ _ C?7» + gi C??igi - 1 

Cl122 " 1122_ c^-^-c^^ + i' 

12 2Z i (C77 i -^)(C^i + l) 12 ' 

45 - - Jg.A'-'fA> e 



Since the potential is transparent, the fermion-breather transmission amplitude is a pure phase factor which can 

(i) (i) 

be read off from c 1122 or cZ 1122 . Interestingly, the result is independent of Qi, 

T (i) = 0]i + gj tViZj - 1 / 4 o^ 
CrH-ZitmZi + V [ ' 

and hence coincides with the result for a baryon with the same values of Zi,rji. 

We now turn to the bound state spinors for the two discrete states which we shall label with double superscripts, 
(the j-th bound state of breather i). In the case of the DHN baryon, the easiest way to get the bound state 
spinors is to compute the residu of the continuum spinor at the poles of the transmission amplitude in the complex 
(■-plane. This works here as well (the poles are at the same location), but the resulting states are not orthogonal and 
therefore not directly suited for the TDHF calculation. The difference comes about because in the breather case, the 
discrete states are not eigenstates of the Hamiltonian, so that there is an ambiguity which linear combination one 
should choose, as also noticed by DHN Q. We impose two criteria: We demand that the bound state spinors tp^'^ 
satisfy the orthogonality and normalization conditions 

dxipV'M ( x ,t)i>^(x,t) = S jk , (44) 

and that the (local) scalar condensate matrix is diagonal, 

^)(x,i)^ (i ' fe) 0M) = 5 ]k 4> {h3] {x,t)^\x,t). (45) 



If one insists on Eqs. (|44| I45[) . one finds that the spinors are unique up to overall phases and interchanging the state 
labels (1,2). This choice then enables us to use the TDHF equation, Eq. @, in the standard form for the discrete 
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states as well. Our results agree with DHN who have reasoned somewhat differently in the semi-classical path integral 
approach. The resulting bound state spinors can be represented in the form 

MM = + /a W) v W + /W>(y«)V« + /^'V^^) 2 . (46) 

The factor % in the lower spinor component accounts for the different transformation properties of left-handed and 
right-handed spinors under Lorentz boosts, leading to more symmetric coefficients. These coefficients are found to be 





. i(*,2) _ 
— J2 — 


e 122 — 


J 112 


= z,-(i - 


ft), 


e 2 


= / (4 ' 2) = 


J*- 2 ) - 
e 112 — 


J 122 


= l + Q l 




e 112 


= '4 (m) 


— e 2 


_ f(i,2) 
— J 122 


= -(1- 


ft), 


e 122 


— J2 




f(i>2) 

— J 112 


= Z S ;(1 + 


0,). 



(47) 

The normalization factor can be chosen as real and positive, 

w / *ft(^ 2 -i) y /2 f481 

All the results of this section reduce to the corresponding single baryon results if we set Qi = — 1. The (real) product 
Vi^V® can be identified with {/W of the static baryon in that case. The spinors agree up to irrelevant overall phase 
factors, ipv* 1 ' goes over into the negative, if>( % ' 2 ^ into the positive energy bound state spinor. Both in the baryon and 
the breather case, ip( hl > and ^^ ,2 ' are related by charge conjugation which reads 

^(i,2) = 75 ^(i,i)* (49) 

in our Dirac basis. 

Consider the computation of the fermion density next. Assuming the lower bound state to be fully occupied and 
an occupation fraction v for the upper bound state, the total fermion density of breather i is given by 

p = i/J\ty (M)t V (M) . ( 5 °) 

The fermion density of the lower bound state is cancelled exactly against the fermion density induced in the Dirac 
sea. If one simply inserts the bound state spinors (I46[) into expression (|50[) . one gets a rather intransparent result. A 
more convenient way of accessing the fermion density is as follows. Since the vector current for each single particle 
state is conserved, it can be represented in terms of a pseudoscalar field P as 

d^ = — ► r = ^ v d v P. (51) 

We use the convention 

(0 = f? "iV M=f_? nV (52) 



We note in passing that the axial current can then be expressed as follows, 

jg = e^j v = d»P, (53) 

however it is not conserved in the GN model. Clearly, P is only defined up to an additive constant. By a judicious 
choice of this constant, one gets a very simple expression for P of breather i, 

p( i) = (Vi V2 1 (54) 

The best way of computing the fermion density and the fermion current is to insert this expression into (1511) . i.e., 

p« = d x P^, jW = -<9 t P w . (55) 
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The normalization condition then reduces to 

dxp {l) = P (l \x = 00)- P (i) {x = -00) = 1. (56) 



Indeed PW of Eq. (|54[) has the shape of a kink evolving from —1/2 to 1/2 with our choice of the additive constant. 

Finally, let us recall how the self-consistency condition ([6]) arises, dropping the breather label for the moment and 
returning to the original parameters (e, b). The scalar density from the continuum spinors yields 

V'cV'c = (^eV'f)i + bPc^ih 

= -^7Y S ( 57 ) 

Upon using the vacuum gap equation, the first term gives self-consistency. The 2nd term can be summed over the 
Dirac sea (integration over <i£ with the appropriate measure) with the result 

d( f C 2 + ± \ (j 1 \ ^ / 2 \ 2V 1 V 2 

\rw ) ( ^ c)2 = "(TT^y y 1 « arctane J -D- (58) 

The discrete spinors contribute 

= T ^ = 2 _Y^1 (59) 



to the condensate. Denoting the fermion numbers of states 1,2 by N,uN, respectively, the condition that (|47|) and 
(|48| cancel reproduces exactly the DHN result, Eq. ((6]). 

The breather is a type II TDHF solution, like the DHN baryon, in the classification scheme of Ref. |6j . Incidentally, 
one can easily check that the breather is in general not a self-consistent solution for the massive GN model. In the 
massive model only the baryon limit Q = — 1 is self-consistent [16(. Hence we do not expect the DHN breather to 
survive in the nonrelativistic limit. It is a genuine relativistic object, like the kink or antikink. 



III. BREATHER-BREATHER SCATTERING 

Here we follow almost literally the strategy proven successful in previous studies of baryon-baryon scattering Q . 
We write down an ansatz for the two-breather problem, starting from the known single breather solution. Next, 
"reducible" parameters are determined from the asymptotics of the incoming and outgoing breathers. The remaining 
"irreducible" parameters can be found by solving the Dirac equation algebraically. This procedure yields a reflec- 
tionless, time dependent scalar potential, together with the continuum and bound state spinors. We then check the 
self-consistency of the solution and compute the fermion density. The final results will be illustrated with a few 
examples. Technically, the present problem is more involved than baryon-baryon scattering which it generalizes to 
scattering of collectively excited hadrons. We believe that the problem is nevertheless worth the effort. When going 
from baryon-baryon scattering to scattering of any number of composite bound states, we observed a kind of factor- 
ization which allowed us to reduce the dynamical multi-baryon problem to baryon-baryon scattering [9] . Similarly, we 
expect the solution of the breather-breather scattering problem to be the key element in a future study of scattering 
of any number of breathers, baryons and composites thereof. 



A. Ansatz and asymptotic conditions 

We are now dealing with two (boosted) breathers with unavoidable notational complications. In Sect. UH the basis 

exponentials referring to breather i were denoted by V^, see Eq. ([37]). The same exponentials will be used in the 
scattering problem for breathers 1 and 2. However it turns out that both the notation and the bookkeeping are 
somewhat simpler if we label the exponentials of breather 1 by Vi,2 and those of breather 2 by V3 4 from now on, 

Vl = exp { 2z[ ~ ^7 J ' V2 = Vl ' 

^3=exp(^-^) I V^Vi. (60) 
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We shall always assume v\ > v 2 in the following. The scalar potential is written as 

S = ^, (61) 

where numerator and denominator are multivariate polynomials of all four exponentials Vk- The degree of these 
polynomials is simply determined by multiplying the scalar potentials of two distinct breathers and keeping all the 
terms appearing there, but with unknown coefficients. Thus the basic structure of Af and T> consists of 25 terms 
generated from 

{i + v 1 2 + v 1 v 2 + vi + vivi)(i + vi + v 3 v 4 + vi + viv'i). (62) 

The ansatz for the numerator Af then reads 

Af = 1 + an Vi + a 22 Vi + a 33 Vf + a 44 Vi + a 12 V 4 V 2 + a 34 V 3 V 4 

+a 1122 vivi + a 1133 vivi + a 11M vivi + a 2233 vivi + a 2244 vivi + a 3344 vivi 

+ai 233 V 4 V 2 Vi + a 1244 V 4 V 2 Vi + a 1134 ViV 3 V 4 + a 2234 ViV 3 V 4 + a 1234 ViV 2 V 3 V 4 
+ 0112233 v(vivi + 0112244 vtvivt + ^113344 V(V£Vl + a 223344 viv 3 'vi 

+a 112234 ViV 2 2 V 3 V 4 + a 123344 V 4 V 2 ViVi + a 11223344 ViViViVi . (63) 

The denominator T> has the same structure as Af with all a-coefficients replaced by 6-coefficients. Similarly, the 
continuum spinors are parametrized as 

« P (g ~ I 



Here, Afi,Af 2 have the same polynomial form as Af, Eq. (|6"3"|) . with all a-coefficients replaced by c- and d-coefficients, 
respectively. Like in baryon-baryon scattering we expect four distinct bound states, two for each breather. They will 
be labelled by 1 and 2 (for breather 1) and 3 and 4 (for breather 2), where the identification of the bound state with a 
particular breather always refers to the asymptotic region. The polynomial structure for bound states asymptotically 
belonging to breather 1 can be obtained by multiplying the numerator of the bound state of breather 1 with the 
numerator (or denominator) of S of breather 2 (20 terms), 

(Vi + v 2 )(i + v x v 2 ){\ + ^ + v 3 v 4 + vi + vivf). (65) 

This suggests the following ansatz for the first bound state spinor, 

,(i) ( Af^\ 

Here, 

A/f) = e [ 1) V 1 + e [ 2 ] V 2 + e { H 2 ViV 2 + e&ViVf + e^ViV? + e^V^Vi 
+e[ 1 4 ) 4 V 1 Vi + eg^Vi + eg\v 2 V 3 V 4 + e^V^ 2 
+e[\ ] 233 ViV 2 V 3 2 + e { H 234 V 2 V 2 V 3 V 4 + e{\{ 44 ViV 2 Vi + e^^ViV 2 
+e[ 1 2 ) 234 V 1 V 2 V 3 V 4 + eiW^ViVi + e^V^Vi + e^ i4 V 2 V 2 Vi 
+c ( iL^ViV 2 V 2 Vi + e ( i 2 \ 334i V 4 V 2 V 2 Vi, (67) 

and Af 2 nas the same structure with all e^-coefficients replaced by /(^-coefficients. Similarly, bound states related 
to breather 2 have the structure generated by a bound state of breather 2 and the numerator (or denominator) of S 
of breather 1, 

(V 3 + V 4 )(l + V 3 V 4 )(1 + Vi + V ± V 2 + V 2 2 + ViVi). (68) 
The 3rd bound state spinor then becomes 
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with 



A/} 



(3) 



4 3 V 3 

, (3) 
T e 223"2 



J3), 



v?v 3 



+eil223^2 2 ^ 

+e<§ 33A V x V 2 V£V A 

J3) 



„(3) T/ 2 T/ , (3) v V 2 
- e 334 l/ 3 V 4 + 63441/31/4 

A2 T r2i 



e 11224^1 2 ^2 2 ^4 



e 11334^1 2l/ 3 2 ^4 



(3) 
e 12344 



TOW? + e 2 3 2 ) 3 34 V 2 2 V-3 2 V4 + 



22344^2 V3V4 



+ e 1122334^ / l 2 ^ / 2 2 ^ / 3 2 ^4 + &\\ 223 ± 4 V l V 2 V 3 V4 



(70) 



Once again has the same structure with all e^-coemcients replaced by /(^-coefficients. The two remaining 

bound state spinors can simply be obtained by charge conjugation, 



(71) 



Our labelling is such that in the baryon limit, bound states 1, 3 become the negative, 2, 4 the positive energy discrete 
states. This ansatz leaves us altogether with a large number of coefficients to be determined, namely 2x25 for S, 
2x25 for the continuum spinors and another 2x20 for each bound state spinor. Even if we make use of Eq. (jTTj) we 
need to determine as many as 180 coefficients. 

The next step consists in reducing this large number by exploiting the asymptotic conditions in the initial and final 
states where the breathers are well separated. This is again analogous to what has been done for baryons in Refs. [1,0]. 
Consider S first. Initial and final breathers can be projected out by letting (Vi, V 2 ) or (V3, V4) go simultaneously to 
or infinity. The relationship between S(Vx, V 2 , V3, V4) for breather-breather scattering and S(Vi, V 2 ), S(V 3 , V4) for the 
two individual breathers is 

lim S(yi,V 2 ,V 3 ,Vt) = S{V U V 2 ), 

V3,V4— >o 

lim S(Vi,V 2 ,V 3 ,V 4 ) = S(a 1 V 1 ,a 2 V 2 ), 

Vsj.Vi— yoo 

lim S(y u V 2 ,V a ,V4) = S(a 3 V 3 ,a 4 V 4 ), 

Vi,V2-»-0 

lim S{V x ,V 2 ,V 3 ,Va) = S{V 3 ,Vi). (72) 

Vl ,V2— >oo 

As illustrated in Fig. [51 the (complex) factors <Ji account for the fact that during the collision the breather undergoes 

S(a 3 V 3 , 04 Vi) SViVl.c^vy 

v u v 2 - 




V 3 , Vi -> Q V u Vt -> 00 



S(Vi,V 2 ) S(V s ,Vi) 

FIG. 6: Asymptotic reduction of scalar potential for breather-breather scattering. Time runs in the vertical, space in the 
horizontal direction. Outgoing breathers 1 and 2 experience a time delay and a shift in their phase, as indicated by the complex 
factors Oi. 



a time delay and a possible change in phase. Thus they encode the whole asymptotic scattering information. Since 
S is real, they have to satisfy 

02 = a* , 0-4 = 03. (73) 
In the baryon limit, U\ — V\V 2l U 2 = V3V4 and the Oi are related to the (real) time delay factors 5\ 2 introduced in 

t— = cno-2 = |o"i| 2 , -!- = 0-3CT4 = |o- 3 | 2 . (74) 

012 021 
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FIG. 7: Like Fig. [(J but for continuum spinors. The incoming spinor 2 and the outgoing spinor 1 acquire transmission amplitudes 
r' 1 ' 2 ' from scattering on the other breather, since they are incident from the left in this figure. 



The condition <5i2<$2i = 1 found for baryons has been confirmed here for breathers in the form 

Cl<T2<73<T4 = 1- ( 75 ) 

The asymptotic reduction for the continuum spinors tp^ is similar to the one for except that one has to take 
into account that the spinor acquires a transmission amplitude in some cases, see Fig. [7] Since the scattering wave 
functions are defined to be incident from the left, the transmission amplitudes affect the two breathers on the right 
side of Fig. [71 namely incoming breather 2 and outgoing breather 1. Denoting the transmission amplitude of the 
spinor due to breather i by Tw as in Sect. Ill D| we get 

T^ c {a 1 Vi,a 2 V 2 ), 

T^^(V 3 ,V 4 ). (76) 

Unlike the a. k which are determined while solving the Dirac equation, the T^ 1 ' are already known from the single 
breather problem, see Eq. (j4*3"|). 

Inspection of the ansatz (|63l) shows that these asymptotic relations relate all coefficients to single breather coefficients 
and oVs except for those multiplying V?Vf, V?V?, V%Vg, V^VZ , V x V 2 Vi, V^Vf , V?V 3 V 4 , V r 2 2 V r 3 V 4 , ViV 2 V z V A . 
This reduces the number of parameters in S and tpc, from 2x25 to 2x9. This reduction is less effective than for 
baryons, where the corresponding reduction is from 2x9 to 2x1 parameters, but nonetheless useful for solving the 
Dirac equation algebraically. These various numbers for baryons and breathers have actually a simple interpretation, 
depicted in Fig. [8l Since our ansatz ()63[) has been derived by multiplying 2 polynomials for single breathers, the 
coefficients fit naturally into a square matrix. From Eq. (|76[) . we see that the first and last row and the first and last 
column of this matrix are filled with reducible parameters, whereas the innermost 3x3 square contains all irreducible 
coefficients. For baryons, the corresponding square is only of size 3x3 and the irreducible inner part a single element. 
This explains the above counting of reducible and irreducible coefficients. 

B. Results for scalar potential and continuum spinors 

As explained in detail in Refs. @, Q, the unknown coefficients can be determined algebraically by inserting the 
ansatz for S and the spinors ipQ into the Dirac equation. Here, we immediately turn to the results for breather- 
breather scattering. We first report on the complex scattering factors Oi describing how Vi in one breather is affected 
by the collision with the other breather. General properties of the <Ji have been given in Eqs. (|73ti75|) . We find that 
the <Ji factorize as follows, 



r lim il>dVi,V2,V 3 ,V4.) = 

3 jV4— >0 

lim ^ c (V ll V2,V3,V 4 ) = 

lim J C (V 1 ,V 2 ,V 3 ,V 4 ) = 

lim ipdVuV^VaM) = 



01 — 013#14, °2 — #23^24, 
03 = #31#32, 04 = #41#42- 



(77) 
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Vi VWi vi 




1 




Vi,V 2 -4 




v? 








VM 


V 3 , V, 
-> 


irreducible 


V 3 , Vi 

— > QC 


vi 












Vi , Vi — > oo 





FIG. 8: Square matrix into which the coefficients of A/", £>, A/i , Mi for breather-breather scattering can be naturally fitted. The 
innermost 3x3 square contains the irreducible coefficients, all other coefficients are asymptotically reducible as indicated in the 
figure. 



The individual factors 9ij may be interpreted as scattering amplitude of V, on Vj , where each Vi stands for a kink or 
antikink constituent of one of the breathers, and have the values 



713 



'14 



'23 — 



7 24 



1 _ ViZ 2 - V2Z1 

1 _ m + t] 2 z 1 z 2 

O41 r\\-ri2Z\Z 2 

1 _ 7?1^1^2 + ??2 

6*32 r\\Z\Z 2 - r/2 ' 
1 _ i]iZi - rj 2 Z 2 
6*42 r/iZi +r) 2 Z 2 ' 



(78) 



Just like in the baryon case, one can actually understand the results (|77l78p in simple terms. First notice that the 
bound states can be identified via the singularities of the fermion transmission amplitudes, i.e., the poles of T^ 1 -* and 
T^ 2 \ Eq. P5|) . in the complex (-plane located at 

Z X 1 

si — — ) C2 — -- 



Vi Z]_r\Y 

(3 = —, Ci = (79) 
r\2 Z 2 r\ 2 

Expressing the 8ij through these complex (-values yields the simple expression 

sj t <,i 



This enables us to identify the scattering factors <7; with transmission amplitudes evaluated at complex £ values, i.e., 
for bound states. First we rewrite using the variables (j rather than 77.;, so that the pole structure is manifest, 

C + C1C + C2 



C-CiC-C: 



2 



T(2)(c) = c+ac+a (81) 

S — S3 S — S4 
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A comparison with ((77)) then shows that 



01 

02 
03 
0"4 



(r^ccx))- 1 , 

(T (2) (C 2 ))^, 

T (1) (C 3 ), 

T (1) (C 4 ). 



(82) 



Thus, like for baryons, one can relate the two-body scattering data to single breather input, namely the transmission 
amplitude for bound states, evaluated at complex spectral parameters. Notice also that the total fermion transmission 
amplitude for the two-breather system factorizes into T^T^ with the symmetric result 



no = n 



(83) 



Since all of these asymptotic quantities are independent of Qi, they are actually identical in the breather and baryon 
cases (for the same Zi,i}i). 

Equipped with the scattering factors 6ij, we can now present the results for all the coefficients in a concise form. 
Although the asymptotic conditions yield only the reducible coefficients along the periphery of the square shown in 
Fig. [8j it turns out that all coefficients except the center element (i = j = 3) can be generated by a simple algorithm 
from single breather input. Thus only the coefficients 01234, &1234, C1234, c/1234 are truly irreducible, exactly like in the 
baryon case. Note that the reduction we are now talking about became only apparent after solving the problem and 
was not anticipated by us, not being related to asymptotics. Nevertheless, it is very useful for presenting the results 
for a large number of coefficients in a compact fashion. 

Let us first explain how to generate all the coefficients except the ones multiplying V1V2V3V4. Using the matrix 
scheme shown in Fig. [8j we introduce the 5x5 square matrix 



M = 



Z 1 
1 

1 

1 

Vi 



n2 n2 A n n n q2 n2 n2 n2 r\ 
t/ 31 t/ 32 P3lf 32*741 P42 ^4X^42 t/ 31 t/ 32 tf 



q2 r{2 
7 31 17 32 

#32 
^31^32 

1 



732^42 


931^41 
1 



2 (\2 
41< 7 42 

#42 

#41 #42 

#ll 
1 



/ 31 L/ 32 l/ 41 l/ 42 \ 
(\2 (Y2 
17 32 t7 42 

931#32#41#42 
9 2 fl2 



J Z1 U 41 
1 



(84) 



Consider the numerator M of S first, and write down the numerators AA 1 ) and J\f^ of S for the individual breathers 
1 and 2, cf. Eq. (|38|) . Convert the polynomials A/^ 1,2 -* into vectors according to their monomial content, using the 
same basis as in Fig. [51 



/ 1 



ul 1 ) 



a (1 V 2 



\ 



ViV 2 



22 v 2 



( 1 

(2)i 



ui 2 > 



V 2 

11 v 3 



a 

a$V 3 V* 



" 2 2 v 4 



(85) 



All coefficients appearing here are single breather coefficients given in Eq. (|39|) . Then, multiplying the matrix M with 
U^ 1 ' from the le: 
coefficient 01234, 



U^ 1 ' from the left and ui 2 -* from the right yields Af for breather-breather scattering, correct except for the single 



Af = Ui 1} MUi 2) . 



(86) 



The one missing, non-trivial coefficient 01234 will be given explicitly below. 

Somewhat miraculously, the same procedure works for the ^-coefficients of the denominator T> of S, as well as for 
the c, d-coefficients of the numerators A/1.2 of the continuum spinors. The matrix M is always the same, so that 
this representation of the results is very economical indeed. In all cases, the vectors U can be constructed from 
the single breather results in a fashion analogous to what was done in Eq. (|85|) . This boils down to replacing all 
the aW coefficients by the corresponding &W,cW,dW, see Eqs. (|39|42p . The resulting vectors will be denoted by 
, Uc , accordingly, and we get 



V 
Ai 
A" 2 



U'^MU^, 
U^MU^. 



(87) 
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These relations yield the correct coefficients except for 61234, C1234, 1^1234. 

This completes the presentation of all the coefficients entering the scalar potential and the continuum spinors for 
breather-breather scattering, except for those multiplying V\ V2 V3V4. The 4 missing coefficients CI1234, 61234, C1234, c?i234 
are the most complicated ones, not related to single breather input in any obvious way. Actually, we have found 
that they are proportional to the irreducible baryon-baryon coefficients af 2 , bf 2 , cf 2 , df 2 of Rcf. [8] with a simple 
proportionality factor depending only on the Qi, 



2Qi 2Q 2 
ai234 ~ Q\ + IQI + I a ^ 



This equation remains valid if we replace a by b, c, d on both sides, with exactly the same proportionality factor. The 
irreducible baryon coefficients have been spelled out in Ref. @ and are somewhat lengthy. To keep the present paper 
self-contained, we list them here once again, using the more structured form given in Ref. Q, 



B 

a 12 = 


a\a\ 
, (P1P2 - 
di2 


VB12), 


bf 2 = 


b\b\ B 
—. — i>i2, 
d\2 




c 12 — 


c l c l / 
«12 




df 2 = 


012 


-B 12 ), 



(89) 



with one-baryon coefficients 



2(Zj + 1) 

Zi(z? + iy 

Zf + V 

(zf + ixcm-Zi^mZi + iy 

2(2z?-e v ?(zt + r)) 
(zf + i^cvi-z^ViZi + iy 

zf-i 

Zf + V 

Zf-l 



Pi = 



Zt + l-2ZfC V f^ 

(z? i)vK 2 



and two-baryon coefficients 

dii = 



1 2zf-evf(zt+iy 

2(rj 1 Z 1 - r) 2 Z 2 )(r) 1 Z 2 - rj 2 Zi){rj 1 ZiZ2 + r? 2 )(r? 2 Z X Z 2 + rji) 
^1(^-1)^-1) 

2zfzUv i + v 4 2)-vW2(zf + i) (zl + 1) 

V W 2 (Z$-1)(Z$-1) 



(90) 



12 — ^2^2/74 1U74 i\ • 



This close relation between breather-breather and baryon-baryon scattering was unexpected, and we do not have a 
simple explanation for the proportionality (|88p . 

Summarizing this subsection, there seems to be no major additional complication when going from baryon-baryon 
to breather-breather scattering, in spite of the increased number of parameters and the more complicated dynamical 
processes which are now described. Our strategy of solving the Dirac equation via the ansatz method works equally 
well in both cases. When we tried to present our results for the large number of coefficients in the breather case 
in the most compact fashion, we discovered surprisingly close connections between breather-breather scattering and 
individual breathers (for reducible coefficients), but also between breathers and baryons (for irreducible coefficients). 
These insights may give us hints as how to simplify and possibly generalize this whole calculation in the future, 
perhaps by exploiting more efficiently the integrability of the GN model. 
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C. Results for bound state spinors 



The bound state spinors can be found either by ansatz and the solution of the Dirac equation, or by analytic 
continuation from the continuum spinors. In both cases, one has to take linear combinations which satisfy the 
orthogonality condition and the additional conditions on the scalar density discussed for the single breather bound 
states, Eqs. (I44I45|) . There are four bound states, two for each breather. The relevant ansatz has already been given 
in Sec. MI Al We will again try to present the e, /-coefficients for states 1 and 3 in the most efficient way, looking for 
factorization of our results. The states 2 and 4 can then be obtained for free, using charge conjugation, Eq. ([7Tj) . 

To present the results for the spinor of bound state 1 in a compact way, we introduce the following vectors, 





1 - c 






w 2 * = 



\±(1 + Q 1 )V 2 2 / 
Then the numerators entering Eq. (|66j) can be represented as 

A/; (1) = Zi^W+MUf 
A/" 2 (1) = FiW^MUf 



/ \ 
1 



l + Qi 

V ±(1 - Q^V? J 



(92) 



C=Ci 



C=Ci 



- Z 4 V 2 W+M t£ 



(2) 



C=C2 



C=C2 



(93) 



The normalization factor for the spinor of bound state 1 can be inferred from the single breather at times where the 
two breathers are well separated, 



(94) 



with Cg 1 ^ from Eq. 



Similarly, the results for the spinor of bound state 3 require the vectors 



/ ° \ 

i-q 2 



\±(1 + Q 2 )V 4 2 J 
The numerators in Eq. (|69p are then found to be 

M 3) = e 13 e 23 z 2 v 3 w+M T 14 1 

13 6 23 V 3 W 3 M T 







1 + Q2 



(95) 



V ±(1 - Q 2 )Vi ) 



9u0 24 ViW;M T Ul 1 ) 



C=C3 



M 3) 



C = C3 



014^24^2 V 4 W+M T U% 



(1) 



C = C4 
C = C4 



The normalization factor for the spinors of bound state 3 is 



C< 8 > = 



(96) 



(97) 



with Cq 2 ^ from Eq. 



The extra factors 9ij in (|9^|) as compared to (|9"3"]) must be a consequence of the definition 



80]) which distinguishes between i > j and i < j. 



D. Self-consistency and fermion density 



The ansatz method yields a time-dependent, transparent scalar potential for the Dirac equation, with the correct 
boundary conditions for breather-breather scattering. We still have to verify self-consistency so as to be sure that we 
have found a TDHF solution of the GN model. The way self-consistency arises always follows the same pattern. The 
scalar density for the continuum spinors can be decomposed into 



2C 



C 2 + i 



-S. 



(98) 



19 



The first term gives self-consistency as usual, whereas the 2nd term can be represented as linear combination of the 
scalar densities of bound states, 

(A^ch = a (1) ^ (1 V (1) + a (! W (3) - (99) 

Upon integrating over d£ and making use of the self-consistency conditions for the individual breathers with occupation 
fractions i/ 1 ' 2 ) , we find analytically 



2ir V 2C 2 



2ir \ 2C 2 



a™ = -(l-i/ x >), 

= -(I-*/ 2 )). (100) 



This is equal and opposite to the contributions from the discrete states and therefore cancelled exactly if we sum 
over continuum and bound states. Like baryon-baryon scattering, breather-breather scattering is a type III TDHF 
solution. 

Consider the fermion density next. We expect that the situation is the same as for baryon-baryon scattering, i.e., 
the total, vacuum subtracted fermion density can be expressed in terms of the bound state densities p^ = and 
p<® = p^ as 

Ptot = AT(>yi) + l ,( 2 y 3) )- (ioi) 

The total fermion number is 

Nf = N (v^ + z/ 2 )) . (102) 

Like for the single breather, the density for a bound state can most conveniently be represented in terms of a 
pseudoscalar field P, see Eqs. (|51ti55l) . The result for bound state 1 in breather-breather scattering has the following 
simple form, 

P (D = v ,(Dt^(i) = JLZl 

H v v dxV 1 

Ti = iu[ 1) diag(-l,0,0,0,l)MUi 2) . (103) 

Similarly, we find for bound state 3 

T 3 = iu«Mdiag(-l,0,0,0,l)u£ 2) . (104) 

Identity (jlOll) is then best verified numerically. This close relation between the total fermion density and the density 
from the bound states is the reason why we discussed the bound states in some detail. Eqs. (|10Hll04j) will be needed 
in the following to compute the fermion density in breather-breather scattering analytically. 



E. Illustrative results 



We have presented the full, analytical solution of breather-breather scattering in the GN model, exact in the large N 
limit. The asymptotic scattering information displays complete factorization, and can in fact be predicted on the basis 
of the single breather. This is a manifestation of the integrability of the GN model. The (reflectionless) fermion-double 
breather transmission amplitude is simply the product of two independent fermion-breather amplitudes, Eq. (I83|) . It 
has poles in the complex plane of the spectral parameter £ corresponding to four bound states, two per breather. The 
asymptotic breather-breather scattering information is encoded in the complex numbers <Zj which multiply the basis 
exponentials Vi in the outgoing channel, see Figs. [B] and [7] As usual for solitons, they lead to a time delay, but also to 
a change in the phase of the breather oscillation. They can be inferred by evaluating fermion-breather transmission 
amplitudes for one breather at the complex spectral parameters (fc corresponding to a bound state of the other 
breather, see Sect. UlLBl This generalizes similar findings for baryon-baryon scattering in Ref. [8[ to breather-breather 
scattering. 
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The TDHF solution does not only yield the S-matrix, but also the full time evolution of the fermion wave functions 
in the whole space. Since the mean field is a classical quantity, we can illustrate breather-breather scattering in 
the same way as in classical, nonlinear systems like the sine-Gordon model. There is nothing wrong in prescribing 
simultaneously the initial position and velocity of the breathers, as they are not subjected to the uncertainty principle. 
The novel aspect of these dynamically generated breathers as compared to classical bosonic theories are the fermions. 
Fermions have been treated in a fully quantum mechanical manner, using the TDHF approach. They populate the 
breathers, and their dynamics can be followed in detail by monitoring the fermion density. Although the polarization 
of the Dirac sea is fully taken into account, we have seen that the total fermion density can be computed using only 
bound state spinors and occupation fractions. Since the self-consistent TDHF potential in the GN model is always 
reflectionless, one expects that the kinks and antikinks making up the breather will repel each other, whereas the 
fermions should always move forwards. This behavior has already been observed in the case of baryon scattering. 

Since this paper is somewhat technical, let us first summarize how a full scattering calculation can be done with 
CA systems like Maple or Mathematica. This should help the reader to navigate through our rather complicated 
set of equations. One first has to set up the matrix M, Eq. (j84|) . the vectors U^ 2 ^ d , Eq. (|85|) and the vectors 
Wf, Eqs. (|92I95|) . The required coefficients are the single breather coefficients which can be found in (|39[) for a, b- 
coefficients and in (|42l) for c, d-coefficients. Eqs. (I86I87I93I96P can then be used to generate conveniently the basic 
polynomials Af, V entering the TDHF potential S, A/1,2 of the continuum spinors and N^ 1 ^ of the bound state spinors 
for bound states 1 and 3. This construction misses the irreducible terms which have to be added by hand, using the 
a, b, c, d-coefficients from (|88ti91[) . The normalization factors for the bound states are given in (|48l94l97p . Together 
with the basis exponentials Vk, Eq. (jHO)) . these ingredients enable us to construct the mean field S from Eqs. (|6HI63[) . 
the continuum spinor ip^ from the bound state spinor ?/> (1) from (12511271) and V (3) from (JBHTD]). The fermion 

density can be readily evaluated using Eqs. (|10Hll04j) . This procedure may sound a bit complicated, but we remind 
the reader that we have to communicate 180 coefficients, each consisting of several factors. We have not found any 
simpler ways of making our results available. 

One can actually apply these results to various related problems. Choosing the two breather velocities differently, 
one deals with the breather-breather scattering problem in an arbitrary Lorentz frame. Choosing the same velocity 
also yields a solution, a (marginal) bound state of two breathers with mass equal to the sum of the two constituents. 
The distance between the breathers can be arbitrarily chosen, like for baryons and their composites [3, 8]. If we set one 
of the Q-parameters equal to —1, the corresponding breather reduces to a baryon and we can handle breather-baryon 
scattering and bound states. If we then tune the Z parameter of the baryon so that it becomes a widely separated 
kink-antikink pair (Z — > i or e — > 00), we can send one of the kinks to spatial infinity and account for breather-kink 
scattering and bound states. This encompasses in particular the analog of the breather-kink bound state that has 
been called a "wobble" in sine-Gordon theory [UGil, a topologically non-trivial breather. Finally, if we set both Q's 
equal to —1 we recover baryon-baryon, baryon- kink and kink-kink dynamics. This is nothing new but can serve as a 
useful test of the algebra. 

Since the results are in closed anlaytical form and their implementation is straightforward using CA, there is no 
point in showing many concrete examples here. The best way of illustrating these time dependent results is anyway 
by means of animations of the scalar potential and the femion density. This is easy on a computer using Maple or 
Mathematica, but unfortunately impossible here. We therefore only discuss briefly two examples of breather-breather 
scattering for illustrative purposes. 

In the first example we choose the breather with parameters e = 2, A = 2.8 as incident from the left and the 
breather with e = 2.0, A = 1.1 of Fig. [5] as incident from the right with equal and opposite velocities (771 = r/^ 1 = 2). 
Fig. [5] shows the time evolution of the scalar potential. One recognizes the oscillations of the incoming and outgoing 
breathers as well as the collision region. Fig. [TU] gives a view at the fermion densities of these colliding breathers. It 
is hard to understand what exactly is happening from these plots. We therefore give a few snapshots of the scalar 
potential and the fermion density during the collision in Fig. 11. The fat curves are S, the thin solid curves are the 
fermion density from bound state 1 and the dashed thin curve the fermion density from bound state 3. Since these 
are separately conserved, it makes sense to split up the total fermion density in this manner. This enables us to 
keep track of the individual lumps of fermions which indeed cross each other, each one moving only in the forward 
direction. 

In our 2nd example we collide two breathers with a more pronounced kink-antikink structure with parameters 
e = 700, A = Vl + e 2 and e = 700, A = 400. The boost parameter 771 = I/772 = 1.005 now corresponds to a low 
relative velocity so as to be able to illustrate the "breathing" and the interaction dynamics on a similar time scale. 
Fig. [T2] shows the scalar potential, Fig. Q2] the total fermion density during the collision. Here, complicated things 
happen with the fermions which seem to disappear and reappear elsewhere in an almost discontinuous fashion. To 
further explore what is going on during the collision, in Fig. [14] we present again a sequence of snapshots. In the 
first and last frame, the two fermion bound states are attached to their "home" breather, as expected asymptotically. 
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FIG. 9: Time evolution of scattering of breathers with parameters e = 2, A = 2.8, v = 0.1735 (breather l)and e = 2, A = 1.1, v = 
0.6753 (breather 2), during collision with rji = I/772 = 2. 




FIG. 10: Time dependence of fermion density for the collision process of Fig. [9] 



Inbetween however, we see that the fermions hop from one kink or antikink to the next, spending most of the time 
near the nodes of the potential. Kinks and antikinks approach each other up to some minimal distance where they 
are reflected. At this point of closest approach the fermions tunnel through, as can be seen in an animation by 
the disappearance of one peak in the density and the simultaneous reappearance of another peak elsewhere. These 
tunneling processes are behind the apparently discontinuous behavior of the density in Fig. 1131 In the collision region, 
the two individual bound state densities of states 1 and 3 do not belong to any particular breather. In the 3rd frame 
of Fig. 1141 bound state 1 is located at the outer kink-antikink pair, bound state 2 at the inner kink-antikink pair. 
This is the (temporal) midpoint of the collision. In the neighboring frames the fermions are attached to 2 or even 
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3 different kinks and antikinks. This figure allows us to follow the hopping mechanism from kink to kink in detail. 
Similar hopping processes were also observed in kink-antikink scattering processes 0, @] ■ Clearly, they are intimately 
connected to the fact that the self-consistent potential is transparent. 




FIG. 12: Time evolution of scattering of breathers with parameters e = 700, A = \/l + e 2 , v — 0.3634 (breather 1) and 
e = 700, A = 400, v = 0.6362 (breather 2), during collision with rft = 1/772 = 1.005. 
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IV. SUMMARY AND OUTLOOK 



In this paper, we have extended previous works on baryon-baryon scattering to breather-breather scattering in the 
massless GN model. We have first introduced the breather originally due to DHN and illustrated its behavior in 
space and time with a few examples. It can be regarded as a vibrational excited state of the baryon of semi-classical 
type, characteristic for the large N limit. We then had to recast it into a form better suited for the scattering 
problem. We used an ansatz method developed in the context of baryon scattering problems. The basis exponentials 
entering the joint ansatz for scalar potential and spinors in the TDHF approach have to be complexified, otherwise 
everything goes through as before. In this way one can derive exact, analytical expression for breather-breather 
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scattering with arbitrary initial conditions and parameters. This contains also breather-baryon or baryon-baryon 
scattering as special cases, as well as bound state problems if the velocities of both scatterers are chosen to be equal. 
The calculation involves a large number of coefficients which have to be determined partly asymptotically (from 
single breather input), partly by solving the Dirac equation algebraically. We have made every effort to present the 
results in the most compact form. While doing this, we encountered a number of simplifications and observed more 
algebraic structure than anticipated. We do not yet fully understand these features, but they point to the possibility 
of simplifying further our calculation. Right now, somewhat ironically, the computations with Maple are perhaps the 
lesser problem as compared to the task of presenting the results in a digestable form. 

Since breathers are kind of exotic objects, at least in particle physics, one may ask whether it is worth the effort 
to study them in such great detail. We think that if one is interested in solving the GN model as completely as 
possible, there is no way around considering breathers. In this context, it is instructive to look at the simpler problem 
of the classical sine-Gordon equation for a moment. There, one may ask what is the most general multisoliton 
solution corresponding asymptotically to spacelike separated individual solitons (i.e., disregarding solutions with a 
finite density of solitons like soliton crystals). This has been answered by inverse scattering theory some time ago. The 
result is a known algebraic solution consisting of any number of interacting solitons and breathers, all with different 
velocities [ijj]. If we ignore the breathers for a moment, this can be compared to the general multi-baryon solution 
of the GN model discussed in Ref. Q . The main difference is that in the GN model, two baryons can have the same 
velocity if they have different fermion numbers, thus describing multikink-bound states absent in the sine-Gordon 
model. Since breathers also arise in the GN model, it is plausible that the most general TDHF solution of the GN 
model (with a finite number of solitons) will also consist of breathers and solitons. In the baryon problem, we found 
that the only input needed to tackle the n-baryon problem are the solutions of the one- and two-baryon problems. 
This was interpreted as signature of factorization on the level of the composite states. Likewise, one might expect that 
a similar generalization exists with breathers, and that the two-breather problem solved here is sufficient to deal with 
any number of breathers, using again factorization. If correct, this would indicate that the two-breather scattering 
problem solved here may actually play a central role and take us a long way towards the most general TDHF solution. 
As a side benefit, we would also be able to generate the most general transparent potential of the Dirac equation, a 
problem similar to the one which has been solved for the Schrodinger equation long ago by Kay and Moses [201 ] . 
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